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| ' We have introduced coherent state neutral pion into Antisymmetrized Molecular Dy- 

namics. With the aid of coherent state technique, it becomes possible to calculate transition 
matrix elements of the pion field operator and to study excited states containing pions. For 
' large pion-nucleon coupling f n N > 1.6, pions have a finite expectation value and bring large 

y—( ' energy gain in 12 C. We discuss two aspects of pionic effects in spectroscopy; the LS inter- 

action like effect and the mixing of different nucleon parity states, which would modify low 
£SJ ■ energy nuclear levels. 

> : 

in . 

OO ■ §1- Introduction 

O ; 

■ Nuclei have been basically understood as nucleon many-body systems in which 

nucleons move in a mean field and interact via small residual interactions. In shell 
models, each nucleon single particle wave function is assumed to have its own orbital 
angular momentum, spin and parity (/j7r), provided that the spherically symmetric 
mean field consists of central and spin-orbit (LS) parts. This basic picture of nuclei 
has been successful in describing low-lying states of most of medium to heavy nuclei, 
while there are several exceptions such as clustering states in light nuclei. 

On the other hand, the long range part (r > 2 fm) of the bare nucleon-nucleon 
. £^ potential is described by the one-pion exchange potential (OPEP) having strong 

^ ■ tensor part, which mixes different partial waves. This mixing is essential for the 

binding of deuteron, but makes it difficult to treat exactly in many-body systems. 
Thus tensor force has been usually treated in the form of effective central and LS 
forces in solving nuclear many-body problems in spite of its importance in bare 
nuclear force, hoping that the explicit role of the tensor force is not large in nuclei. 
Actually, the first order effects of tensor force vanish in the Fermi gas state, due to 
the cancellation in (2 J + l)-weighted sum. 

If pions have expectation values in nuclei or in nuclear matter, the cancellation 
does not work and tensor force may play a dominant role. There have been a 
lot of discussions on the possibilities of pion condensation in nuclear matter. For 
some time, it was considered that strong nucleon-Z\ short range repulsion might 
suppress pion condensation to emerge, provided that the Landau-Migdal parameter 
follows the "universality", g' NA = g' NN ~ 0.6, and that it is density-independent P 
Recent observations of the non-quenching in Gamow- Teller giant resonance sum 
rule^ clarified that the short range repulsion between nucleons and A resonances is 
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not very strong, g' N ^ < 0.25,^ suggesting that the pion condensation would come in 
exist at least in high-density nuclear matter. Also in recent ab-initio calculations of 
light nuclei (A < 12) with realistic bare AW potentials, it has been shown that the 
OPEP contribution is dominant in the total potential energy.™ This result suggests 
that the cancellation in the Fermi gas is not working well in actual nuclei and that 
it would be necessary to consider the explicit role of pions more seriously. 

On these backgrounds, pion condensation and explicit role of tensor force in 
nuclei has now attracting renewal interests. In a relativistic mean field framework, 
it is demonstrated that neutral pions can condensate in the surface region of nuclei, 
and that this condensation enhances the binding energy of jj closed nuclei, such 
as 12 CP In their model, single particle states are first prepared to have fixed Ijn, 
and those states having the same j are mixed to gain the potential energy from 
pions. This mixture also plays a role of LS like potential!^*™ Since pions mix 
different parity (but same j) states, the yrast single particle states (the lowest energy 
single particle states for a given j, s 1/2 ,p 3 /2, d 5/2 , / 7 / 2 , S9/2, ^ii/2> H3/2) wil1 nave tlie 
largest energy gain. On the other hand, those states having smaller j at around 
the Fermi energy will be pushed up from the mixing with the lower energy single 
particle states. It is interesting to note that the last four yrast single particle states 
are in charge of the nuclear magic numbers of 28, 50, 82 and 126. More recently, 
Charge and Parity Projected Hartree-Fock (CPPHF) method has been developed in 
order to take account of the coupling of proton and neutron single particle states 
generated by OPEPP^ It has been shown that the charge projection enhances the 
tensor contribution by around three times in the case of 4 He nuclei. 

At this stage, it would be desirable to extend the scope of pion and tensor 
force study from the ground state and single particle states to nuclear excited level 
spectroscopy with specified J w , which has richer information on wave functions. 
Specifically, we are more interested in constructing a framework in which explicit pi- 
onic degrees of freedom are incorporated, rather than introducing tensor interaction, 
since we believe that it is more fundamental to describe nuclear many-body system 
with pions. 

In this work, we introduce coherent state pions^ into Antisymmetrized Molec- 
ular Dynamics (AMD/^JOJ} anc i discuss the pionic effects on excited states of 12 C. 
Pion coherent state enables us to calculate matrix elements of the pion operator 
with different states as well as the expectation value with a given state. In AMD, 
nuclear wave function is represented by the Slater determinant of nucleon Gaussian 
wave packets, which is wide enough to describe clustering states as well as shell 
model states. By using the product of nucleon AMD state and pion coherent state, 
we can evaluate the transition matrix element of the Hamiltonian operator contain- 
ing nucleon and pion operators. Therefore, it becomes possible to project the wave 
function to the eigen state of given J n and to diagonalize the Hamiltonian matrix 
consisting of wave functions with different nucleon and pion configurations. 
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§2. AMD with Coherent State Pion 



The nucleon-pion basis state is assumed to be the product of the nucleon AMD 
stated E3 and the pionic coherent state}*' 

\9(Z, /)) = |1Pamd(Z)) |<M/)> • (2-1) 
AMD wave function is a Slater determinant of nucleon Gaussian wave packets, 
|^amd(Z)) = AH\i> Zi MXi) (Z = { Zi ;i = l,2,... A}) , (2-2) 



2v \ 3/4 

— exp [-i/(r - .z/v^) 2 + z 2 /2] 

IT J 



(2-3) 



where IxfxD represents spin-isospin wave function. 

Pion coherent state introduced by Amado et alP is represented as, 



|<M/)) = exp 



dkf{k)a\k) 



|0) 



(2-4) 



By setting the commutation relation of the annihilation and creation operator, a{k) 
and a)(k'), as [a(k),a^(k')] = 5(k — k'), we can show that the above pion coherent 
state is an eigen state of the positive frequency operator 



he dk 



-.d{k)e 



ft + \r,t) , 

J y/(2^2u k /c 

ft+Xr,t)\$M)=<Pf(r> t )\*M) 



tpf(r,t) 



he dk 



f(k)e 



ikr—iLo^t 



(2-5) 

(2-6) 
(2-7) 



Since the bra state is an eigen state of the negative frequency operator, </>(~>(r,t), 
and the pion operator is a sum of and c^ - ** , we can easily calculate the transition 
matrix element of the pion operator as follows, 



(r,t) = ^ + \r,t) + ^-\r,t) , ft~\r,t) 



r,t) 



(2-f 



($M)\kr,t)\$A9))=MAf,9) x fcf(r,t) + <p g (r,t)) , 



j dkf(k)g(k) 



(2-9) 



(2-10) 



Now we consider the following Hamiltonian of TV-body nucleon and pion system 
containing the second quantized pion operator <f> in the axial vector P-wave pion- 
nucleon coupling, 



N o 
Pi 



1=1 i<j 



1 

2hc 



dr 



V0(r) • V(/)(r) + ^f(r 
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N 

+ Ttifo-VO&fi) , (2-11) 

where fi n = m^c/h, and we have omitted the time dependence in the pion part. 
The matrix element of this Hamiltonian is evaluated as 

n = I^zjwW^ = Hn{2 ' z ' ] + nAL 9) + nM2 ' z '' 1 9) ' (2 ' 12) 

KM, 9) = J ^[{V(^(r) + ^(r))} 2 + ^^ /(r) + 999(r) }2] > (2 . 13) 

■H 7rN (Z,Z'J~,g) = ^ [ drS(r)-V& f (r) + <p g (r)) , (2-14) 

, , = QPamd(S)| ZiViTQiSir-n) \Vamd(Z')} 
[ ' " (*amd(^)|^amd(Z')) ' 1 ; 

Here Wat is the usual AMD Hamiltonian matrix element including iV./V interaction. 

In the actual calculation, we have expanded the pion eigen function <p(r) in 
local Gaussians, whose centers and amplitudes are the variation parameters. Thus 
we can apply the cooling equations for these pion parameters and nucleon phase space 
parameters Zj's. We have made the non-relativistic approximation in the calculation 
of pionic state norm, 



MrCfjS) - exp 



dr<ff(r)ip g { 



(2-16) 



The Hamiltonian form of Eq. (|2T1[) with pion-nucleon P-wave interaction is 
the simplest one. In addition to coupling with charged pions, higher dimension 
terms such as pion-nucleon S- wave interaction coming from N<p 2 N coupling would 
have visible contributions when pions have large expectation values. Charged pions 
should give similar energy gains to neutral pions, but coherent state treatment of 
charged pions mixes different charge states, and this charge mixing may lead to seri- 
ous problems in the discussion of excited levels. In order to overcome this problem, 
it is necessary to perform the coupled channel calculation of different nucleon and 
pion charge states or to perform isospin projectionP2) These are beyond the scope 
of this paper, and will be discussed elsewhere. On the other hand, higher dimension 
terms such as the pion-nucleon S'-wave interaction are not expected to give large 
contributions in energy, since the number of pions is around 0.1 in 12 C nuclei in the 
present framework as is shown later. 

§3. An Example of Application — 12 C Nucleus — 

We have applied AMD with coherent state pion to the study of 12 C nuclei. In 
the ground state of 12 C, nucleons occupy the single particle states of 0si/ 2 and 0p 3 / 2 , 
both of which are the yrast single particle states, then the pionic effects are expected 
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to be large. While the 3a cluster model generally describes the excited levels of 
this nucleus very well, the first excited state (2+) is calculated to be too low. Since 
the "spin-orbit" splitting of the 0p 3 / 2 and 0pi/ 2 single particle states is responsible 
to the 0^-2+ level spacing, pionic effect to push up the 0pi/ 2 level may appear as 
the increase of E*(2f) in 12 C. It would be also interesting to study un-natural 
parity levels such as 0~, 1 + , 2~, . . ., whose excitation energy might decrease due to 
the coupling to the natural parity nucleon state with 0~ pionic state. 

In this paper, we apply the simplest model of AMD with coherent state pions 
as the first step. We show the results with projection after variation (PAV); we 
first construct the intrinsic state by using the cooling variational method for the 
parametrized wave function of Eq. (|2-lj) . and projection to specified J n has been 
carried out from the prepared intrinsic state. We find that the effects of parity 
projection before variation in 12 C nuclei are not large when pions are included ex- 
plicitly, while it has been found to be important for spectroscopic studies of light 
nuclei without pions PE) i n the cooling stage of intrinsic energy, the imaginary part 
of <ff(r) is a redundant degree of freedom as is clear from Eqs. ()2-13|) and (|2-14|) with 
(ff = ip g , then we cannot control the imaginary part which is given randomly in the 
initial state of variation. Thus we have limited the pion eigen function y?(r) to be 
real. We find that there are many local minima especially for small f n N values, then 
we have selected the lowest energy wave functions from several candidates obtained 
from different random seeds. Since we do not take account of the isospin projection 
which enhances the tensor force effect by around three times, we use larger f n N value 
in the range f^N = 1 ~ \/3- For finite nuclei, we use the scaled coupling constant 
f n j^^ = y/(A — l)/Af w N m order to include approximately the effects of the Fock 
(exchange) term of OPEP, which requires quantum corrections in a field description. 

As the effective nucleon-nucleon interaction, we start from Brink-Boeker type 
two range Gaussian interactions, which approximately reproduces the binding ener- 
gies of 4 He, 16 (28.4 and 128.9 MeV, BBOl^J) or the binding energy of 4 He and 
nuclear matter saturation (28.1 MeV and E/A = -16.8 MeV at p = 0.165fm~ 3 , 
BB02). Since the potential energy from pions is very large, it would be important 
for nucleon-nucleon interaction to have saturation property in order to avoid col- 
lapsing. The interaction range is chosen to be shorter than that of, for example, the 
Volkov interaction. When we include pions, we employ the nucleon-nucleon inter- 
action BB07T, which is a little modified from BB02 to reproduce the ground state 
energy and the first excited state energy of 12 C. Parameters of these interactions are 
summarized in Tabled 

Table I. Volkov and Brink-Boeker-Okabe interaction parameters. Central interactions are given as 
the form, V = £- = i v t exp(-r 2 /M?)(l ~ M i ~ M,P CT P T ). 





Mi (mi) 


vi (MeV) 


Mi 


[Xi (fm) 


v 2 (MeV) 


M 2 


/ttJV 


Volkov 


1.6 


-83.34 


0.575 


0.82 


144.86 


0.575 





BBOP 1 


1.2 


-253.798 


0.2186 


0.6 


924.631 


-1.551 





BB02 


1.2 


-258.3 


0.25 


0.6 


950.00 


-1.658 





BBOtt 


1.2 


-256.0 


0.25 


0.6 


950.00 


-1.658 


1.63 
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Fig. 2. Composition of the pionic energies for f^N = 1.63. Thick solid, long-dashed, and short- 
dashed lines show total, nucleonic, and pionic energies, respectively. Among the pionic energy, 
kinetic and interaction energies are also shown by thin solid lines. Filled circles show the 
minimum points of Etot, E w , and E(0 + ). 



In Fig. we show the intrinsic state energy as a function of the pion number 
expectation value, 



«M/)|«M/)> * — J dr ^ r) ^ r) ■ <M 
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At small firN values around one, pure nucleon state is energetically favored. When 
we increase f n N, pion- nucleon interaction gives very large binding, and the optimal 
state has finite pions for f^N > 1.6. The total pionic energy amounts to be around 
—90 MeV in the case f^N = 1-63 as shown in Fig. [2 In well developed pionic states, 
the nucleus loses energy in the nucleon part TLnn instead of gaining pion-nucleon 
interaction energy efficiently. This feature is similar to the case of pion condensation 
in high density nuclear matter.^* 




Fig. 3. Pion number dependence of the total energy (top), energy difference from 0^ state (middle), 
and the nucleonic abnormal parity probability (bottom) in 12 C. Filled circles show the energy 
minimum points for each J* '. 



In the upper panel of Fig. [3 we show the results of total energy after J n projec- 
tion from the cooled intrinsic wave functions under n n constraint. All the J 71 " states 
have their minima at finite n n when we adopt f n N = 1-63. It is interesting to find 
that natural parity states favor smaller n n , and un-natural parity states favor larger 
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Finite number pion is expected to act as the LS-like interaction and to increase 
the excitation energy of 2f . In the middle panel of Fig. |21 we show the energy 
difference E(J W ) — E(0j) as functions of n w for f n N = 1-63. At = where the 
present model is equivalent to the normal AMD, 2f has small excitation energies, 
which is a feature of a cluster models. At around n w ~ 0.05, l\ state starts to 
go up. This increase of energy difference at > 0.1 is not a consequence of the 
nuclear shrinkage, but the result of pionic LS-like effect. Actually we find that the 
calculated rms radius grows in the region n w > 0.07. 

Contrary to the positive parity rotational states, 0j~ and if states go down as 
the pion number increases. This is due to the coupling to the nucleonic abnormal 
parity states such as, 

|#(0-)) = IM0-)) + IM0 + )> ® I<ZV((T)) , (3-2) 

for the 0~ state. In order to demonstrate this point, we show the nucleonic abnormal 
parity probability P^° n ' in the bottom panel of Fig. 13 At zero pion number n n = 0, 
all the states should be purely described in nucleonic state (P^ bn ' = 0), but the 
probability goes up to around 5 % and 3 % at the projected energy minima for 0~ 
and 1 + states, respectively. Other rotational levels are also calculated to contain 
the abnormal nucleonic parity component of around 1 %. If these are true, it may 
be interesting to observe pion knock-out reaction, which leaves the nucleus in the 
un-natural parity states. 

Table II. Energy components in 12 C levels. For + state, calculated total energies are shown in 
the parentheses. All the energies are shown in MeV. 





r 




E* 


T N 


Vc 


Vcoul 


Vls 




Volkov 


0+ 


( 


-92.4) 


234.3 


-320.6 


8.9 


-15.0 






2+ 




4.2 


234.5 


-320.6 


8.8 


-11.1 






4+ 




12.0 


235.4 


-317.5 


8.8 


-7.2 






3" 




19.5 


244.5 


-314.1 


8.7 


-12.0 




BBOtt 


0+ 


( 


-92.4) 


187.5 


-228.8 


7.2 


0.01 


-58.2 




2+ 




4.6 


190.9 


-233.3 


7.2 


-0.03 


-52.5 




4+ 




10.6 


196.3 


-224.9 


7.1 


-0.15 


-60.1 




1+ 




17.0 


198.1 


-195.3 


7.0 


0.25 


-89.6 




tn 




27.0 


204.0 


-183.5 


7.0 


0.48 


-89.1 




3" 




8.7 


194.7 


-228.9 


7.2 


-0.04 


-56.5 



The ground state and the first 2 + state energy can be reproduced in AMD 
without pion effects when we adopt strong LS interactions,^ but the wave functions 
in these two descriptions are very different. In Table [Hj we compare the energy 
components for 12 C levels in AMD with Volkov interaction with strong LS interaction 
(Vls = 1800 MeV) and in the present model with moderate LS interaction (Vls = 
900 MeV). The energy difference of + and 2 + mainly comes from the LS interaction 
in the case without pions, while the pionic energy is the main source of the energy 
difference when pions are included. In addition, it is interesting to find that the 
LS interaction acts in the reverse way — LS interaction is more attractive for + 
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without pions, but it weakly acts repulsively for + with pions. 

§4. Summary 

In this paper, we have developed a new framework to include pionic degrees of 
freedom in the nuclear many-body systems, Antisymmetrized Molecular Dynamics 
(AMD) with coherent state pion. Compared to the mean field treatment of pions, 
the present model has a merit that we can evaluate the transition matrix elements 
of the pion-nucleon coupling term containing the second quantized pion operator 
4>. This enables us to calculate the excitation energies for specified J n with explicit 
pion degrees of freedom, through the parity and angular momentum projection from 
the intrinsic state. It is also interesting to find that the pion coherent state has a 
norm, and the pionic state overlap (<P 7T (f)\<P- K (g)) reduces the total state overlap, 
(¥(Z,f)\¥(Z',g)) : 

We have applied this model to the study of 12 C structure. The LS'-like effects 
of pions can be seen as the increase of the 2f state excitation energy. It is also 
suggested that explicit pionic state |<?tv(0 + )) <g) |$„.(0 - )) can admix to the 0~ state 
with around 5 % when we adopt f w N = 1-63. 

In order to obtain firm conclusions on the explicit pionic effects in nuclear 
structure, further theoretical and experimental developments are mandatory. First, 
charged pions should be included in the framework. Combined with the charge 
(isospin) projection, it is expected to enhance the tensor effect by three times. We 
have simulated this enhancement by increasing the pion-nucleon coupling f^N, but 
the exchange of proton and neutron may lead to non-trivial effects which cannot be 
mimicked by increasing f w ^. Next, we should take care of the exchange term (Fock 
term) and the zero range (5 type) part of the OPEP in a more reliable manner than 
to scale the coupling constant by y/(A — 1)/A times. Including the Landau-Migdal 
interaction may be an efficient way for this problem. Extension of the wave function 
is also an important direction. Since we have assumed that the total wave function 
is a product of nucleonic state and pionic coherent state, nucleonic part of the wave 
function is common to the zero and one pion states. Thus the nucleonic part of the 
wave function has to contain both of the T = and T = 1 states, which may lead to 
the underestimate of the binding energy. Works in these directions are in progress. 
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